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In a previous paper we suggested that, for a given p mode, the excitation function is the same as the component of the solar 
background noise that has an identical surface spherical harmonic projection (over the corresponding range of temporal 
frequency). An important consequence of this surmise is that the excitation of overtones of a given angular degree and 
azimuthal order will be correlated in time. In this note, we introduce the basic principles and a mathematical description 
of correlated mode excitation. We use simple, illustrative examples, involving two modes. Our treatment suggests that 
in the real observations, any signatures of the correlation would not appear as a correlation of the output amplitudes of 
overtones, but rather as subtle modifications to the power spectral density at frequencies between the central frequencies 
of the overtones. These modifications give a contribution to the observed peak asymmetries. 

© 0000 WILEY- VCH Vcrlag GmbH & Co. KGaA, Wdnhcim 



(N 



Oh 
I 

O 



> 

oo 
m 
m 
m 

^' 

o 

oo 

o 



X 



1 Introduction 

It has now been well established that solar p-mode peaks 
observed in the frequency power spectrum have profiles that 
are asymmetric in frequency. The modes are excited stochas- 
tically in a thin layer at top of the convection zone, and the 
asymmetry is believed to result from important characteris- 
tics of this excitation (e.g., see Nigam & Kosovichev 1998; 
Nigam et al. 1998; Kumar & Basu 1999a, b; Severino et al. 
2001; Jefferies et al. 2003). First, there is a contribution to 
the asymmetry from the spatial location, radial extent, and 
multipole properties of the acoustic source. Second, there is 
a contribution from what is commonly referred to as cor- 
related background noise. This background noise is signal 
from the convective granulation. Since the convection gives 
rise to the acoustic source that excites the modes, the granu- 
lation and modes are correlated. These correlations give rise 
to asymmetry. It is the observational consequences of such 
correlations that are the subject of this paper 

In Toutain, Elsworth & Chaplin (2006), we presented 
a framework to explain the contribution of the correlated 
noise to the asymmetry. We hypothesized that the excitation 
function of a mode of angular degree /, azimuthal degree m, 
and frequency i/q, is the same as that component of the so- 
lar background (granulation) noise that has the same spheri- 
cal harmonic projection, ¥;,„, in the corresponding range in 
frequency in the Fourier domain. It is this component of the 
noise background that gives a contribution to the asymme- 
try. 

An important implication of the above is that overtones 
of a given (l,m) should have excitation functions that are 
correlated in time. Indeed, correlation of the excitation fol- 
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lows naturally from invoking correlations with the back- 
ground noise. CoiTelations are not expected between dif- 
ferent components of the same mode. That is because the 
Yi„i for {l,rn) and (l,—m) are orthogonal, and are therefore 
assumed to have independent, i.e., uncorrected, excitation. 

Our main aim in this paper is to introduce the basic prin- 
ciples and a mathematical description of correlated p-mode 
excitation. The paper is didactic in nature, and uses sim- 
plified examples - involving just two modes - to illustrate 
the principles. The paper should be viewed as providing the 
necessary background, and framework, needed to consider 
the more complicated case of the full spectrum of overtones 
of a Sun-like oscillator. We leave a description of how the 
effects are manifested in the solar p-mode spectrum to an 
upcoming paper. 

It is important to stress that correlation of the excitation 
in time does not imply correlation of the mode amplitudes in 
time. This can be understood by considering the analogy of 
damped, stochastically forced oscillators. Modes of differ- 
ent frequencies will be 'kicked' by the common excitation 
at different phases in their oscillation cycles, and provided 
the frequency differences are at least several times greater 
than the peak linewidths, there will be significant differ- 
ences in how the amplitudes vary in time due to the excita- 
tion, a point we return to in Section|2]below. The correlation 
is instead manifested in a more subtle manner in the fre- 
quency power spectrum: specifically, it modifies the power 
spectral density of the spectrum at frequencies between the 
central frequencies of the overtones, where the wings of the 
Lorentzian-like mode profiles interact, and gives a contri- 
bution to the observed peak asymmetry. It is tempting to 
say this another way: that "the background is modified" 
between the modes. However, that would be a misleading 
statement, in that with correlated background and correlated 
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excitation the clear distinction between modes and back- 
ground no longer exists. The next generation of 'peak bag- 
ging' mode-fitting codes will need to account for this blurred 
distinction. 

The layout of the rest of the paper is as follows. We be- 
gin in Section|2]by showing that if two modes are excited by 
the same functions in time, this does not always imply the 
oscillation amplitudes will also be correlated in time. We 
then turn to the issue of the mode profile asymmetry, begin- 
ning in Section |3] with a recap of the discussion in Toutain, 
Elsworth & Chaplin (2006) on how correlated background 
noise gives rise to mode peak asymmetry. This is followed 
by the main part of our paper: In Section 2] we show how 
correlation between the excitation of different modes modi- 
fies the observed power spectral density (analytical descrip- 
tions are provided in AppendixlAJi. We finish in Section|5]by 
discussing the implications of these modifications, flagging 
work in progress. 

2 Time correlation of amplitudes of different 
modes 

In order to illustrate the impact of correlated mode excita- 
tion on correlations in time between the output amplitudes 
of different modes, we made simulations of stochastic har- 
monic oscillator timeseries. The Laplace transform solution 
of the equation of a forced, damped harmonic oscillator was 
used to generate timeseries of the output velocity of artifi- 
cial modes on a 40-sec cadence, in the manner described by 
Chaplin et al. (1997). The oscillators were excited at each 
time sample with small 'kicks' from timeseries of random 
Gaussian noise. 

For each simulation we generated artificial timeseries 
realizations of two modes. Both modes were excited by the 
same timeseries of kicks, meaning their excitation was 100 % 
correlated in time. The first mode always had a natural fre- 
quency of 3010 /xHz; the frequency of the second mode was 
changed from one simulation to another, and, over the full 
sequence of simulations covered the range 2990 to 3010 /iHz. 
We chose this range to illustrate results for two correlated 
modes, closely spaced in frequency. Both modes always had 
intrinsic damping rates that gave natural linewidths of F = 
1 /iHz for their resonant peaks. This linewidth corresponds 
approximately to the linewidth of modes at the centre of the 
low-l solar p-mode spectrum, and is equivalent to a lifetime, 
r = I/ttF, of ~ 3.7 days. We comment briefly below on 
how the results are affected by the choice of F. 

For the analysis, we used 0.5-day-long timeseries. This 
meant the mode peaks were unresolved in the frequency do- 
main. From each pair of 0.5-day timeseries, we computed 
the RMS velocity amplitudes of the first and second modes. 
We repeated each simulation 100 times, using a new time- 
series of random Gaussian noise to excite both modes. We 
then computed the Pearson correlation coefficient between 
the 100 measured RMS velocity amplitudes of the first and 
second modes. This gave us a single correlation measure for 
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Fig. 1 Correlation of the output amplitudes of two modes, 
as measured from 0.5-day timeseries, as a function of the 
separation in frequency between the two modes. The exci- 
tation of both modes is always 100 % correlated in time. 



a given pair of input frequencies. We then made new sets of 
100 simulations for different input frequencies of the second 
mode. 

A plot of the measured correlations of the output am- 
plitudes is presented in Figure [T] The correlations are plot- 
ted as a function of the separation in frequency, Ai/, be- 
tween the two modes. When the frequency separation was 
zero, both modes had a frequency of 3010 /iHz. When the 
separation was 20 /iHz, the first mode had a frequency of 
3010 /iHz and the second mode had a frequency of 2990 /xHz. 
It is readily apparent that when the frequencies differ by 
more than « 5 /iHz, the Pearson correlation coefficient is 
scattered around zero. Only at small frequency separations 
does the correlation of the output amplitudes become sig- 
nificant (and of course reaches unity for zero separation). 
The rate at which the correlation falls with increasing sepa- 
ration of the input frequencies depends on the ratio Ai'/T. 
Figure [U shows that the correlation coefficient drops to 0.5 
when At^/F ss 2. So an increase of F will result in the 
coefficient of correlation falling off more slowly (and vice 
versa). 

Consecutive overtones of the low-Z solar p modes are ac- 
tually separated in frequency by Av ^ 135 /iHz, which is 
the so-called large frequency spacing. This separation is al- 
most seven times larger than the maximum separation shown 
in Figure[T] If our hypothesis regarding the correlated mode 
excitation is correct - this being that it is the excitation of 
overtones of a given (l,m) that is correlated - the results here 
indicate the measured amplitudes of the modes will not be 
correlated in time. 

3 Recap: impact of correlated background 
noise 

Consider a p mode that is stochastically excited and intrin- 
sically damped. The excitation function is assumed to be 
a random function with zero mean and unit variance (the 
strength of the excitation is folded into the mode amplitude). 
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The excitation function in the frequency domain, E{i'), may 
therefore be written as: 



E^v) = e,.(i/) + iei{u). 



(1) 



Similarly, the noise is assumed to be a random function with 
unit variance, multiplied by a frequency-dependent ampU- 
tude, i.e.. 



^H = \/^KM + ««.M]- 



(2) 



The frequency response of the mode, which we assume may 
be modelled as a forced, damped oscillator of high Q, is 
described as a Lorentzian, i.e., 

where H is the maximum power spectral density (mode 
'height'), X — {h' — i'o)/{T/2), and uq and F are the central 
frequency and linewidth of the mode, respectively. Follow- 
ing, for example, the discussion in Anderson et al. (1990), 
we relate these quantities to determine the solution, in the 
Fourier (frequency) domain, of the oscillator equation de- 
scribing the behaviour of the p mode: 

V{v)^ L{i^)E{i^)+N{iy). (4) 

If we assume that the excitation function and the background 
noise are uncorrelated the limit power spectral density is de- 
scribed by: 

^M = (in^)n = Yf^+"H' (5) 

which includes the usual Lorentzian profile. Here, the an- 
gled brackets ((...)) indicate an average over a large number 
of realizations of the excitation function, and of the back- 
ground noise. 

Let us now assume that the excitation function and the 
background noise of the observed spectrum are correlated, 
so that in each frequency bin 

< CrUr > = < eiUi >= a, (6) 

where a is therefore a coefficient of correlation between 
the background noise and the excitation function. The limit 
power spectral density then takes the more compUcated form: 

H 



P{v) 



1 



1 + 2a^/n/Hx] + n{v) 



(7) 



For non-zero a, this gives an asymmetric profile. 

This scenario is illustrated schematically in Figure |2] 
which we call Case#L The excitation function is repre- 
sented by a timeseries of 'kicks' of random Gaussian noise. 
The kicks excite an oscillator. In the example shown, the 
oscillator has a natural frequency of vq = 2990 yuHz and a 
linewidth of F = 1 /^Hz. When the timeseries of kicks is 
added to the output of the oscillator, and is therefore used as 
background, we will have correlated background noise. A 
direct consequence is that the frequency power spectrum of 
the final timeseries will show a peak with asymmetry. The 
underlying, or limit, frequency power spectrum for Case #1 
is shown in the right-hand panel of Figure |2] 



In the example shown we scaled the timeseries of kicks 
to give a background-to-height ratio of n/ H — 0.5 %. We 
did not add any uncorrelated background to the timeseries, 
so that the background noise was 100% correlated with 
the excitation (i.e., a = 1). Because this correlation was 
positive, the sign of the peak asymmetry was also positive 
(i.e., more power on the high-frequency side of the reso- 
nance). Negative correlation will give a peak showing neg- 
ative asymmetry. 

Finally in this section, we note that for Case#l - and 
for the examples that follow in Section |4] below - the fre- 
quency spectrum of the excitation was white, because the 
timeseries of kicks had a Gaussian distribution in time. This 
is a reasonable, first-order approximation to the solar case, 
since locally the spectrum of the solar granulation is white. 
One could of course be more sophisticated: In the latest 
version of the solarFLAG helioseismology simulator, ar- 
tificial p modes are excited with random noise having a 
more realistic granulation-like frequency spectrum, which 
has power that rises with decreasing frequency (Chaplin et 
al., in preparation; see also Chaplin et al. 2006, for informa- 
tion on solarFLAG). 

4 Impact of correlated excitation of different 
modes 

We have seen that when the excitation function of a mode 
is correlated with the background noise, the profile of the 
mode is asymmetric. Now, we ask what happens to the ob- 
served power spectral density when there is correlation of 
the excitation of different modes. We take the simple case 
of a frequency power spectrum that is comprised of two 
modes. 

As noted previously, consecutive overtones of the low-? 
solar p modes are separated in frequency by ^ 135 /iHz. 
Here, we consider two modes separated in frequency by 
20 /zHz. This smaller frequency spacing exaggerates the im- 
pact of the correlated excitation on the observed power spec- 
tral density, and therefore allows us to show more clearly the 
effect of the correlation. We return again briefly to this point 
in Section |5] below. We only consider the case of 100% 
correlation, again, as an illustrative example of the effects. 
In the real solar p-mode spectrum differences in, for exam- 
ple, the radial dependence of the excitation with frequency 
might be expected to reduce the correlations. The impact of 
changing the magnitude of the correlations is considered, in 
the context of the full solar p-mode spectrum, in a future 
paper (Chaplin et al., in preparation). 

We begin with a simple reference case, which we call 
Case #2. We again assume each mode has a linewidth of 
1 /xHz. The excitation function of each mode is assumed 
to be independent in time. Furthermore, we do not include 
any background noise in the timeseries. The key elements 
of Case #2 are shown schematically in the left-hand panel 
of Figure |3] The timeseries of 'kicks' excite the oscillators 
(modes) giving the outputs of velocity as a function of time. 
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Fig. 2 Case#l: Excitation of a single mode. Correlated background noise is added, which is 100% correlated with the 
excitation function. 



which are then summed to give the final timeseries. The un- 
derlying frequency power spectrum of the final timeseries is 
shown in the right-hand panel of Figure [3] Since the excita- 
tion is independent, the observed frequency power spectrum 
for Case #2 is given by the (incoherent) addition of the fre- 
quency power spectra of the independent oscillator outputs. 

Next, for Case #3, we assume the excitation of the two 
modes is 100% correlated in time. This case is shown in 
the left-hand panel of Figure |4l where the same timeseries 
of random Gaussian noise is used to excite both oscillators. 
Again, we include no background noise in the timeseries. 
The frequency power spectrum of the final timeseries shows 
clearly that the peaks are asymmetric. We stress that this 
asymmetry does not come from the influence of correlated 
background noise (e.g., as in the single-mode case illus- 
trated in Figure|2]i: no background noise was included here. 
Instead, it comes from the interaction of the two modes, 
whose excitation is correlated. An analytical description of 
the power spectral density is presented in Appendix lA.il 

We draw an important conclusion from Case #3: cor- 
related mode excitation can give a contribution to the ob- 
served asymmetry of modes. 

Real frequency power spectra include contributions from 
background. Let us therefore consider what happens when 
we add background noise to the timeseries shown in Case #2 
and Case #3 above. For Case #4 (Figure |5]l, we add some 
independent, uncorrected Gaussian noise to the final time- 
series of Case #2. The addition of this background noise to 
the final timeseries gives rise to a constant background in 
the frequency power spectrum. Since the excitation func- 
tions and the background are all uncorrelated, the observed 
power spectral density is given by the incoherent addition of 
the frequency power spectra of each oscillator in turn, and 
the background noise. 

Finally, for Case #5 (Figure |6l), we add correlated back- 
ground noise. The timeseries of background noise is just 
the timeseries of excitation kicks so that, as in Case#l, the 
background is 100 % correlated with the excitation. We now 



have two factors contributing to the observed asymmetry 
of the peaks. First, there is the interaction between the two 
modes, whose excitation is 100 % correlated. This contribu- 
tion was illustrated in Case #3. Here, we have a second con- 
tribution from the correlated background (like that shown in 
Case#l), which further distorts the mode peaks. An analyt- 
ical description of the power spectral density is presented in 
Appendix lA.2l 

5 Discussion 

Our aim in this paper was to introduce the basic principles 
and a mathematical description of correlated p-mode ex- 
citation. We used simplified examples - involving interac- 
tions between two artificial modes, separated in frequency 
by just 20 ^Hz - to illustrate the principles. We saw that 
when the excitation functions are correlated in time, the in- 
teraction of the wings of the Lorentzian-like peaks gives 
rise to peak asymmetry. We also saw that, provided the fre- 
quency separation of the modes is at least several times the 
peak linewidths, the output amplitudes of the modes will 
not be correlated in time; information on the correlations is 
instead coded in modifications to the shapes of the mode 
peaks. 

For the Sun, we hypothesize that it is the excitation of 
overtones of a given angular degree and azimuthal order 
that will be correlated in time. This follows from the as- 
sumption that the excitation function of a given mode is the 
same as the component of the solar background noise that 
has an identical surface spherical harmonic projection (over 
the corresponding range of temporal frequency). Overtones 
of low-degree solar p modes are separated by ^ 135 /iHz, 
which is somewhat larger than the 20-/iHz spacing we used 
in the simple two-mode examples in this paper When there 
are just two modes, the impact of the correlated excitation 
on the shapes of the resonant peaks will be much smaller 
when the separation between modes is 135 /iHz, as opposed 
to 20 /zHz. However, there are many overtones in the solar 
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Fig. 3 Case #2: The excitation of two modes is independent, and no background noise is added. 
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Fig. 4 Case #3: The excitation of two modes is 100 % correlated, but no background noise is added. 
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Fig. 5 Case #4: The excitation of two modes is independent. Unlike Case #2 the timeseries also includes uncorrelated 
background noise. 
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Fig. 6 Case #5: The excitation of two modes is 100 % correlated, like Case #3; however, here the timeseries also includes 
background noise that is completely correlated with the excitation (see Case #1). 



spectrum, and so even with the larger frequency spacing be- 
tween correlated modes the cumulative effect of the interac- 
tions between all the overtones may be expected to give rise 
to significant modifications to the power spectral density. 
Peak asymmetry arising from these correlations would add 
to the well-known contribution from the correlated back- 
ground noise, and would need to be allowed for explicitly 
in attempts to make inference on the physics of the peak 
asymmetries. We consider the case of the full low-degree 
solar p-mode spectrum in Chaplin et al. (in preparation). 
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A Analytical power spectral density for 
2-mode spectra 

In the following we adopt the naming and symbols conventions 
from Section[3] 



A.l Description for Case #3 

We begin with a description of the power spectral density of Case #3 
(Figure|4}, i.e., a frequency power spectrum comprising two modes, 
whose excitation functions are correlated in time, but there is no 
background noise. Let the frequency response and excitation func- 
tion of the first mode be Li (v) and Ei {v), respectively. The cor- 
responding functions for the second mode are 1/2 (i^) and _E2(i^), 
respectively. The observed complex frequency spectrum is there- 
fore: 



V{i^) = Li{u)Ei{u) + L2{p)E2{u), 
and the observed power spectral density is: 
P{i^) = {\V{v)\'') = [LiMSiM + L2{i^)E2{u)] X 
[Ll{u)El(u) + L*2{u)El(v)]. 



(Al) 



(A2) 



In what follows we shall drop the explicit dependence of the func- 
tions on V. Multiplying out the terms in Equation I A2I above, we 
have: 



P = LiLlEiEl + LiL*2EiE2 

+ L1L2E1E2 + L2L2E2E2 ■ 



(A3) 



Next, we recall from Section [3] that the excitation of a mode is 
formulated in terms of the real and imaginary parts of a complex 
random function. The excitation functions of the two modes may 
be written thus: 



El = eir +ieii, 

and 

E2 = e2r + ie2i. 



(A4) 



(A5) 



We define a coefficient p to describe the correlation of the excita- 
tion of the two modes: 



< eire2r >=< eiie2i >= p. 



(A6) 



The coefficient p is analogous to the coefficient a (Equation |6l(, 
which was defined in Section[3] recall that a describes the correla- 
tion of the excitation of a mode with the background noise. Here, 
p instead fixes the correlation of the excitation of one mode with 
another. 
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The cross term E1E2 in Eauation lA3l mav be written in terms 
of the components of Equations IA4l and lA5l 

E1E2 = [eir + ieii" 



e2r 



te2i 



= eire2r + eiie2i + i[eiie2r — e^ieir] 

= 2p + = 2p. 

The other cross term, E1E2, is also equal to 2p. We may also 
simplify the terms E\E1 and E2E2- Let us take the first of these: 



EiEi = [eir + ieii 
_ 2 , 2 _ 

— Ei^ + Eli - 



[ei,. — leii 
-2. 



2. The expression for P therefore 



We also have that E2E2 
simplifies to: 

P = 2LiLl + 2L2LI + 2p[LiL*2 + LlLi]. (A7) 

We may then expand out the terms (see Section[3} to give: 



P = 



Hi 

1 + xl 



+ 



H2 



+ 2pVlhlh 



1 + X-1X2 



{l + xl){l + xl) 



When p = Q, and the excitation of the two modes is independent, 
we see that Equation I A8I reduces to the sum of two Lorentzians. 
This scenario was described as Case #2 in the paper. When p 7^ 0, 
and the excitation is correlated, the observed power spectral den- 
sity is modified by the third term in Equation I A8 1 which contains 
cross terms between the two modes. The mode peaks then have 
asymmetry. When p = 1, this more complicated scenario corre- 
sponds to Case #3. 



2a 



[ii +12 + 11 + 1*2]. 



(AlO) 



Equation |AT0] above comprises two parts. The first part, which in- 
cludes the terms on the first line of the equation, is just the ob- 
served power spectral density for two correlated modes, i.e.. Equa- 
tions IA7I and IA8I (Case #3) from Appendix I A. II The second part 
has a contribution, n, which is just the power spectral density of 
the background noise, and further terms which describe the corre- 
lation between the modes and the background noise. These corre- 
lated terms modify the observed power spectral density in a non- 
trivial manner, and give a contribution to the observed asymmetry 
of the mode peaks. 

For completeness, we may take one more step and expand out 
the terms in Equation I A 1 01 above to give a final expression for the 
power spectral density: 



P = 



Hi 

l+xl 



(AS) j^^j^ 2py/Hjr2 



l + 2aj-—xi 
til 

1 + 3::ia;2 



H2 



l + 2a 



X2 



{i + xi){i + xiy 



(All) 



When p = 1 and a = 1, Equations I A 1 01 and I A 1 II describe the 
power spectral density of Case #5 (Figure[6ll. 



A.2 Description for Case #5 

Finally, we consider the description of the power spectral density 
of Case #5 (Figure|6j, i.e., a frequency power spectrum comprising 
two modes, whose excitation functions are correlated in time, and 
there is also correlated background noise. The observed complex 
frequency spectrum is now: 

V{v) = Li{v)Ei{v) + L2{v)E2{iy) + N{v). (A9) 

The observed power spectral density is therefore: 

P = 2LiL*i + 2L2L2 + 2p[LiL*2 + LIL2] 

+ NN* + N*[EiLi + E2L2] + N[L*iE*i + LIE2]. 

The first term on the second line of the equation above is just the 
expectation of the noise background described by Equation [2] in 
Section[3] i.e., NN* — n. The other terms in the second line above 
depend on the correlation between the excitation and background 
noise. Let us take the term N* Ei. It may be written in terms of 
the real and imaginary components of A'^* (cf. Equation|2j and Ei 
(Equation I A4b respectively, i.e.. 



iv*£;i = 



n 



- iUiWeir + ieii 



— [eirJir + eiiUi +i{eiinr — eirUi) 



= 2a 



+ = 2q 



The terms N* E2, NE^ and NE2 simplify to give the same ex- 
pression. Putting this all together, the expression for P may there- 
fore be written as: 

P = 2LiL*i + 2L2L2 + 2p[LiL*2 + L*iL2] 
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